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Abstract 

Anisotropic emission of gravitational waves during the merger or formation of black holes can lead to the ejection of 
these black holes from their host galaxies. A recoiled black hole which moves on an almost radial bound orbit outside 
the virial radius of its central galaxy, in the cold dark matter background, reaches its apapsis in a finite time. The 
low value of dark matter velocity dispersion at high redshifts and the black hole velocity near the apapsis passage 
yield a high-density wake around these black holes. Gamma-ray emission can result from the enhancement of dark 
matter annihilation in these wakes. The diffuse high-energy gamma-ray background from the ensemble of such black 
holes in the Hubble volume is also evaluated. 
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1. Introduction 

Accretion onto a black hole (BH) from a non- 
dissipative medium, such dark matter (DM) 
dominated Universe, is rather inefficient, mainly due 
to the absence of a cooling mechanism. However, 
dark matter distribution can be highly modified in 
the presence of a black hole. Enhancement of dark 
matter density and spike formation in the process of 
adiabatic accretion of dark matter onto black holes 
has been studied extensively (Gondolo & Silk 1999). 
Here, we study a different mechanism and study the 
response of dark matter to tlie BHs that are ejected 
from their host galaxies and move in a cold dark mat- 
ter background on bound orbits (Mohayaee, Colin & 
Silk 2008) . Indeed, in this article we review in more 
detail the results already presented in a short letter 
previously (Mohayaee, Cohn & Silk 2008). 

Assuming a zero velocity dispersion for DM and 
a relative motion between the BH and DM, and as- 
suming the flow in the frame of the BH is steady 
and uniform at infinity, the density diverges on the 
downstream symmetry axis. This comes about be- 



cause rings of dark matter, concentric with the sym- 
metry axis, shrink down to points on the axis and 
this loss of dimension manifests itself in a singular- 
ity, i.e. a caustic. This caustic is merely a focal line 
and forms due to axial symmetry and it is not a 
stable catastrophe as classified by Arnol'd (Arnol'd 
1986). Unlike a shock which forms in a dissipative 
medium, the particles cannot be trapped in a caus- 
tic and hence a caustic has very little mass in spite of 
its high density. However, due to the finite primor- 
dial dark matter velocity dispersion, a high density 
wake, rather than a line singularity, forms. In this 
case some of the dark matter particles with veloc- 
ities smaller than escape velocities move on bound 
orbits around the BH instead of arriving from infin- 
ity on hyperbolic trajectories. 

If dark matter consisted of weakly-interacting 
particles, such as those which arise in supersymme- 
try or extra-dimension extensions of the standard 
elcctrowcak model, they would annihilate in pair 
and produce a host of particles among which are 
very high-energy gamma-rays. The flux of these 
gamma-rays depends on local dark matter density. 
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Fig. 1. The trajectory of a test particle in a cold (zero velocity 
dispersion) medium past a perturber, e.g. a black hole. 

Hence, the enhanced density of the wake leads to 
an ever-more enhanced flux of high-energy gamma- 
rays. It is this flux which we aim to determine in 
this work. 

How can black holes attain a velocity relative to 
the DM environment ? And can the velocity disper- 
sion of the background dark matter and the BH ve- 
locity be small enough for the wakes' over-density 
to be significantly high ? 

Black holes (BHs) can gain recoil velocity dur- 
ing their formation if their pregenitor stars collapse 
asymmetrically (Bekenstein 1973) or they can have 
gravitational recoil during their mergers with other 
BHs (Fitchett 1983). The kick velocity can be as 
large as a few hundreds of km/s during formation 
(Bekenstein 1973) and has recently been shown to 
reach a few thousands of km/s for maximally rotat- 
ing BHs, with a particular spin configuration, dur- 
ing the merger phase (Campanelli et al 2007) . 

In spite of the recoil velocity, the density enhance- 
ment in the wake can be extremely low, because the 
velocity dispersion of dark matter is usually very 
high. As dark matter collapses into a halo which 
then evolves by accretion and merger towards a fi- 
nal virialised state, its effective velocity dispersion 
increases. Our galactic halo is assumed to have a ve- 
locity dispersion of a few hundreds of km/s. This can 
be compared with the present primordial velocity 
dispersion of dark matter which is only a few cm/s 
for neutralinos (that of axions is 7 orders of magni- 
tude smaller). However, as we go back in time, dark 
matter becomes less and less clumpy and although 
the primordial velocity dispersion of dark matter in- 
creases linearly with redshift (its present value is 
about 0.03 cm/s for neutralinos), its average veloc- 
ity dispersion in the clumps falls as on average they 
become less massive and smaller. Due to this low ve- 
locity dispersion, many of the BHs with a moderate 
kick velocity can escape from their host haloes at 



high redshifts (Favata et al 2004, Merritt et al 2004, 
Portegies Zwart 2000) . The recoil velocity does not 
depend on the masses of the merging BHs but only 
on the ratio of their masses and on the spin configu- 
rations. Hence, at high redshifts where dark matter 
haloes and escape velocities from them are small, 
only a small recoil velocity is sufficient to set a BH 
on its orbit outside the virial radius of its host halo. 
In addition, at high-redshifts the velocity dispersion 
of dark matter is mostly small and hence the wakes 
could be significantly overdense if the BH also has 
a small velocity. The latter inevitably happens near 
the apapses of the radial orbits. 

The physics of the early Universe is unveiling fast 
(Barkana & Loeb 2001). To explain the quasar popu- 
lation at high redshifts. a large number of early mas- 
sive BHs must have existed (see t.g. Madau & Rees 
2001). These BHs are believed to have formed in 
dark matter haloes and grown by merger with other 
BHs and hence could attain recoil velocities. These 
BHs can travel beyond the virial radii of their host 
haloes and move through colder and colder environ- 
ments. At the apapses passages these BHs come to 
rest and there a substantial density enhancement 
can occur. Dark matter annihilation in these high 
density regions could open a new window into the 
early Universe for high-energy gamma-ray explor- 
ers. 

In Section [2l we review the previous results on 
the calculation of the wake density and obtain the 
density profile in the limit as the velocity of the BH 
goes to zero. In Section 3, we evaluate the time du- 
ration a BH spends around the apapsis of its orbit. 
In Section 4, we evaluate the absolute luminosity of 
a BH in high-energy gamma-ray, as a function of its 
mass and redshift. We also evaluate the boost fac- 
tor due to the BH wake relative to the background 
and also relative to the luminosity of the parent halo 
assuming a NFW profile. In Section 5, we evaluate 
the diffuse background, using Press-Schechter for- 
malism and compare our results with the minimum 
flux from the host haloes. In Section 6, we conclude. 

2. Density of the wake 

We consider the black hole to be a point mass 
moving in a cold and almost homogeneous dark mat- 
ter fluid. Due to this motion a wake forms behind 
the moving object which has a enhanced density 
w.r.t. the background. Indeed, the enhanced density 
of the wake has been suggested as the underlying 
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Fig. 2. The density profile of dark-matter (with vanishing thermal velocity) around a BH of mass M = IOOMq moving with a 
velocity of 1 km/s. The coordinates are shown in Fig. [T]and Fig. |3] 



reason for dynamical friction (Chandrasekhar 1949, 
Kalnaj 1971, Henon 1973). In this section we evalu- 
ate the wake density for the three following cases: (I) 
a medium with zero velocity dispersion, (II) when 
black hole velocity is comparable or smaller than 
the velocity dispersion of its environment and (III) 
when the BH velocity is far larger than the velocity 
dispersion of the dark matter medium. 

2.1. zero velocity dispersion: ctdm = 

In the case of zero velocity dispersion, the density 
of the wake can be evaluated by solving a two-body 
problem (Binney & Tremain 1987). The orbit in po- 
lar coordinates (as shown in Fig. [1]) is given by 
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where R is the impact parameter as shown in Fig. [TJ 
In this case, the density profile can be obtained an- 
alytically, using mass-conservation and can be writ- 
ten in the following convenient form 
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where p is the background density, Mbh is the mass 
of the BH moving with velocity Vbh and the dis- 
tance r and angle 6 are as shown in Figures [T] and 
[21 Evidently a singularity forms on the axis behind 



the BH for 9 ^ tt where the density diverges (see 
Fig. [2]). The density can also be written in terms of 
the cartesian coordinates as, 
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and the profile is shown in Fig.[2l where the coordi- 
nates are shown in Figs. [1] and [31 The situation con- 
sidered here is unrealistic and dark matter particles 
do indeed have a finite primordial velocity disper- 
sion. This case is studied in the next subsection. 

2.2. Density of the wake: (Tdm 

In this case, we assume that the DM particles 
have non-zero temperature and their velocities obey 
Maxwellian distribution. The integral form of the 
density profile can be found by using Jeans' theo- 
rem. The calculation is too detailed to be reviewed 
here (Danby & Camm 1957, Gricst 1988). The fi- 
nal expression for the density enhancement due to a 
moving point mass in a thermal environment is 
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Fig. 3. Contour plot of the density field around a moving BH. 
The plots are made using equation l|10| l for a 10"^ BH moving 
with a velocity of 1 km /s in a dark matter field with velocity 
dispersion of 1 cm/s, and after a suitable normalization for 
graphical presentation, colour coding: red-to-blue represents 
higher density to lower density. 



where 
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a nd p ^ Vbh/o-dm, q = 2G'A/bh/?7'^: 
\J V? -V g-^ (cos^cosA — sinSsinAsini^) , p is the density 
of the environment of the BH, and the distance r 
and angle 6 are as shown in Fig. [T] and Fig. [3l 

2.3. Density oj the wake: Vbh S> ctdm 

When the velocity dispersion of the medium is 
very low in the limit Vbh ^ ctdm , the integral in (O 
becomes highly oscillatory and difficult to evaluate. 
The small velocity dispersion, ctdm , of dark matter 
renders the wake density finite downstream along 
the symmetry axis. For a dark matter velocity given 
by — W = (w — V, V, w), the wake moves to a new 
position, and the density enhancement is given by 
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where cosd 



i^TT^ , tana = — and the angles are 
related using spherical trigonometry relation 

cosV' = — COS0COS0 + sin6'sin0cosa (8) 

and subsequently 



W^2(l-cosV') = 

W'^ (1 + - u)/W cose - w/Wsuv6 



(9) 



Under the condition that Vbh S> (Tdm, one can 
demonstrate, using expression dll) that away from 
the negative symmetry axis (downstream) the den- 
sity approaches that for zero velocity dispersion 
After a suitable re-arrangement of expression ([7]) 
one obtains (see Sweatman & Heggie 2004 and also 
Sikivic & Wick 2002) 



p(MBH,2,r-,6i) = 



2p(^)^B^H 
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where the function /(A/bh, 2, a) is given by 



(10) 
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^ 2GMbh 
X ^1 + coa6 — 2\J\ + COS0 n cosa -f 

We integrate expression (|10p numerically and show 
the density contours in Fig. [3l However, one can 
show analytically that for Vbh ^ ctdm , the density 
pO|) along the symmetry axis [B = tt) attains the 
following maximum value: 
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as long as cr'^^ ^ GMbb/t. That the wake den- 
sity is independent of the velocity of the BH. down- 
stream along the symmetry axis, might be surpris- 
ing. However, we recall that the density is infinite 
there for a zero velocity dispersion and a cut-off to 
this divergence is put by finite dark matter velocity 
dispersion. 

The density enhancement along the symmetry 
axis can be obtained by direct integration of pI7|) . 
The density increases with squared-root of the mass 
of the black hole, grows quadratically with z and 
can be shown to also fall slowly with distance (as 
l/^/r). These results are confirmed by the simple 
expression p^ . 

In the case of zero velocity dispersion, the density 
enhancement given by ^ has a 1/Vbh dependence 
and hence diverges as Vbh — * 0. For a non-vanishing 
velocity dispersion, the density enhancement given 
by the integral (|lGp is independent of velocity along 
the symmetry axis as shown by expression (jl2p (see 
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Fig. d]). Off the symmetry axis the density faUs al- 
most hnearly with increasing velocity, and flattens 
at large velocities as predicted by the exponential 
term in expression (jlOp . 







2.4. Density of the wake : Vbh 



The zone of influence of BH decreases with in- 
creasing its velocity and the velocity dispersion of 
its environment. In the limit as Vbh ~* 0, the den- 
sity profile of the wake ^ reduces to 
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where Erfc is the complementary error function 
and r, is the radius of influence of the BH: r, = 
GMbh/^dm. We emphasis that (fT^ is the limit 
Vbyi ^ of ([n]), and is not a unique density profile 
for stationary BHs. Here, we use (fT3)) only as an ap- 
proximation to ^ for slowly-moving BHs. Fig. [H 
shows the dependence of the density enhancement 
on the BH velocity. When the BH is moving fast 
with respect to the background, a significant den- 
sity enhancement only arises in a small zone around 
the symmetry axis (downstream) of the BH. The 
density enhancement also decreases with increas- 
ing velocity dispersion of dark matter environment. 
The highest density enhancement and largest ra- 
dius of influence are achieved for BHs moving slowly 
(Vbh < cdm) in a cold background. 

Expressions dS]), or ^3]) . have been obtained by 
assuming that the BH velocity remains constant, 
which is questionable in our situation where the BH 
is both slowed down by the pull of its parent halo 
and also by the dynamical friction of the wake it- 
self. Since the BH is most luminous near the apapsis 
where its velocity is very small, we expect our as- 
sumption to be indeed valid near the apapsis. This 
is justified by Fig. [3J once the BH velocity is less 
that the velocity dispersion of the background, the 
constant velocity approximation is valid. 

3. Time spent around the apapsis 

The BH remains bound to its central halo if it is 
ejected with a velocity less than the escape veloc- 
ity (measured from the virial radius). Because it is 
ejected from the centre (and also when with a large 
velocity), the BH is on almost radial orbit. The BH 
initial velocity is set as follows. We assume that the 
halo mass is about 2 x 10^ times the mass of the 



BH (Madau & Rees 2001), keeping in mind that 
the validity of the Magorrian relation (Magorrian 
et al 1998) at high redshifts is yet to be confirmed. 
Thus, for a BH of mass Mbh, the virial radius of the 
halo, from which it was ejected, can be determined 
using Mhaio = 47r/3Avir(2)p(2;)i?vir(^) a-nd noting 
that Avir(z) = (IStt^ + 82a; - 39a;2)/fJ(z) and x = 
n{z)-\&Yidn{z) = fJ„(H-z)V[f^m(l + 2)^ + f^A + 
rifc(l + z)2] (see Bullock et al 2001 for details). Hav- 
ing evaluated the virial radius, we can then evaluate 
the escape velocity from the virial radius of the halo, 

using Voscapo = \/2GMhalo/^vir- 

The ejected BH is slowed down by the gravita- 
tional pull of its parent halo and also by the dynam- 
ical friction of dark matter background as 
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where E = — Vi^/2 -|- GMhaio/^vir is the absolute 
value of the energy with which a bound BH leaves 
the virial radius with velocity Vi. Since dynamical 
friction plays a sub-dominant role in braking the BH, 
the values of ln(A) and the background density p 
marginally affect the value of (fT4|l for a BH in its 
initial outward journey. 

By comparing the dynamical friction force to the 
force of the parent halo in ([14]), we can find the 
range of values of the velocity for which the former 
dominates, give by the inequality 



1^bh(V|h + 2S) < 400GA/bh^ 



(15) 



For a BH with ejection velocity (from virial radius) 
of about half the escape velocity, the inequality be- 
comes 



Vbh < 4 X 10 
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\Me J 
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(16) 



in km/s, for which the dynamical friction dominates 
over the pull of the parent halo. 

Time spent at the apapsis is defined to be the 
time during which the velocity of BH reduces from 
the background DM velocity dispersion to zero (at 
the apapsis), i.e. < Vbh < fDM(-z) where a{z) is 
the velocity dispersion of dark matter in the field 
outside the halo. If the dynamical friction dominates 
over the gravitational pull of the halo, in bringing 
the BH to rest at apapsis. This yields 
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Fig. 4. The density Q for Vbh ^ ""DM a-nd GSj for ^BH S> ctdm shown for a IOOMq BH moving in a cold ( of primordial 
velocity dispersion ctdm = 0.03cm/s) [left panel] and hot (ctjjm = 2km/s) [right panel] environment at high redshift (z=10). 
The right panel also demonstrates the 1 / y/r profile of the density. 



"where Coulomb logarithm is set to unity here and 
is not expected to be significantly greater than this 
value. 

If the pull of the halo is the dominant force in 
bringing the BH to rest then the time during which 
the BH can be considered stationary [hence its ve- 
locity < Fbh < cr-DM{z)] is 
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Ho"wever, the density enhancement is most signif- 
icant when the velocity of the BH and the veloc- 
ity dispersion of DM background are very low and 
hence when the dynamical friction force dominates 
over the gravitational pull of the host halo. Subse- 
quently, A^DF is shorter than Athaio- 

The important issue we have not discussed here 
is that the wake does not form instantaneously and 
the time scale for the formation of the wake has 
to be compared to the time the BH spends at the 
apapsis. This would give us a meaningful estimate 
of the radius of the wake. However, we postpone this 
issue to the forthcoming work and here leave the 
radius of the wake, i?cutoff , as a free parameter. We 
estimate the suitable range of this parameter which 
yields a sufficient luminosity. The criteria here are 



that BH luminosity would dominate over that of 
its central halo and the boost over the background 
luminosity evaluated within the same radius, i?cutoff ) 
be far greater than unity. 

4. Gamma-ray flashes from BHs around 

apapses 

Cold dark matter if composed of neutralinos or 
Kaluza-Klein particles would annihilate in pairs and 
produce a host of secondary products, including en- 
ergetic photons (e.g. see Bertone, Hooper, Silk 2005 
for a recent review). The absolute luminosity^ in 
units of 7S~^ of a BH of mass M at redshift z is: 



L{M, z)= 



in 



r^dr{p{M,z,r)y 



(19) 



where is the Schwarzchild radius, r, is the radius 
of influence of the BH, is the ncutralino mass 
100 Gev/c^), {(jv) is the interaction cross-section 
[which we fix at 2 x 10^^^ cm'^/s] and M-y is the 
number of photons produced per annihilation. Note 
that the integral p9|) is independent of angle for 
stationary BHs. 

We had previously found that the wake den- 
sity of a slowly- moving BH (Vbh < ctdm) is well- 
approximated by the wake density of a stationary 
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Fig. 5. The minimum radius -Rcutoff apapsis required 
for the BH luminosity to dominate over that of its parent 
halo is shown by the lower solid lines for different masses. 
For reference, the minimum virial radius, -Rvir, of the halo 
(corresponding to ejection of a BH of mass IOOMq) is also 
plotted. 

BH. By inserting [keeping only the first term] 
in (|19p we obtain the analytic expression (in units 
of 7 s^^) 



Lbh = 1-3 X 10" 
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for the absolute luminosity of a BH where i?cutoff 
and (TDM are given in the same length units. For 
dark matter with primordial velocity dispersion of 
a — 0.03(1 + z) cm/s, expression (^0)) reduces to 



Lbh = 1.4 x 10^^ (1 + z) 
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for the absolute luminosity of a BH, where i?cutoff is 
in unit of parsecs. We recall that in ([20|) and ([2T|) . 
RcutoS is the radius within which the luminosity of 
the BH is evaluated. 

Next, the BH luminosity is compared to the back- 
ground luminosity, LbGi within the same radius in 
the absence of BH, which is given by 



Lbg = 5.36 X l{r^^{l + zfR, 



cutoff 5 



(22) 



where i?cutoff is in unit of cm. 

The BH luminosity shall also be compared to the 
absolute luminosity of its central dark matter halo 
of mass 2 x lO"' A/bh, assuming it has a NFW density 
profile (Navarro, Frenk & White 1997) 
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where fi(z) nm{i+z)^/[n„,{i+zf+nA+nk{i+ 

z)% 5c = (Avir/3)cV(ln(l + c) - c/(l + c)), and 
c is the concentration parameter, for which we use 
the fit [10/(1 + z)](Mvh./M)-° i3 which agrees weh 
with Bullock et al 2001, Hennawi et al 2007 , and 
is slightly lower than Comerford & Natarajan 2007. 
The halo mass within the virial radius is Afhaio = 

Mvir. 

The absolute luminosity L in unit of 7s~^, 
of a NFW halo of mass Afhaio can then be evaluated 
using ([231) in CH) [after setting M = Afhaio, r. = 
Rv^n:^ f s = 0] and fitted by the following functions 
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A lower limit can be put on the parameter i?cutoff 
in ([2T|) . by requiring that iBH/^Nrw > 1- This is 
shown in Fig. [S] The figure shows that at high red- 
shifts, very small cutoff radii are sufficient for the 
black holes to be more luminous than their host 
haloes. 

However, in comparing the BH luminosity to that 
of its parent halo, we have assumed a primordial ve- 
locity dispersion for DM. This assumption becomes 
less realistic at lower redshifts. The velocity disper- 
sion of dark matter outside haloes has not yet been 
studied in numerical simulations. Observations re- 
port on different values, depending on the environ- 
ment. For example in the local group, the velocity 
dispersion has been reported to be as low as 40 km/ s 
(Karachentsev et al. 2003). Furthermore, it is not 
clear how the velocity dispersion in the field evolves 
with redshift. A high luminosity BH studied here, 
requires a low velocity dispersion environment. Here 
we put bound on this velocity dispersion by requir- 
ing the BH luminosity to dominate over that of the 
background and also over that of its central halo. 



Lbh ^ -^BG and Lbh > ^nfw 



(25) 



using expressions pO|l for the black hole luminosity 
^BH, (|22p for the background luminosity Lbgi and 
([M)) for the luminosity of the parent halo Lnfw- 
The relationship (j25p puts an upper-bound on the 
velocity dispersion. The velocity dispersion of a dark 
matter halo as a function of the redshift and the 
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mass of its orbiting BH (recalling that A/haio = 2 x 
10'* A/bh) is given by 



fhalo- 
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where again n{z) = 0.3(1/(1 + z) + 0.3(1 - 1/(1 + 
z)) + 0.7(1/(1 + zf - 1/(1 + z)))-\ The result 
is plotted in Fig. [6l In this figure, the lower three 
curves are upper-bounds [see expression (|25p ] to the 
velocity dispersion corresponding to BHs of masses 
lOOM©, 10"* Mq and lO^M© as marked on the plot. 
The upper three dashed-dotted (blue) curves are 
plotted for reference and represent the velocity dis- 
persion of the dark matter haloes from which the 
BHs were ejected, given by ([26]) . 

Evidently the maximum velocity dispersion out- 
side the halo would be smaller than that inside the 
halo. Here we have given only an approximate out- 
line and we emphasis that a better method to clarify 
the problem of DM velocity dispersion would be to 
study the velocity dispersion outside DM haloes at 
different redshifts in N-body simulations. Although 
less decisive, analytic studies can also be made 
through for example Press-Schechter formalism. 

Further important constraint on i?cutoff can be 
put by studying the formation time of the wake. 
Indeed, on one hand a finite time is required for a 
wake to form within a certain radius and on the other 
hand the time spent at the apapsis is not infinite. By 
comparing these two time scales, namely time for 
formation of the wake and the time the BH spends 
at the apapsis (during which it can be considered 
stationary) we can better determine the radius of 
the wake. However, we do not expect this radius to 
be much smaller than the minimum values of i?cutoff 
given m FigEl More detailed works on this problem 
is postponed to the forthcoming article and here the 
wake radius is left free parameter. 

5. Diffuse gamma-ray background 

Next, we study in a cosmological scenario, ejected 
BHs near their apapses passages, especially those at 
high redshifts. where the merger rate is higher, the 
escape velocities are lower and the recoiled velocities 
are as large as now. The recoil velocity depends on 
the mass ratio of the BHs and not on the masses 
of the individuals, which indicates that ejected BHs 
are expected to be more abundant at high redshifts. 
Hence the ensemble of recoiled BHs might yield an 
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Fig. 6. The minimum value of the velocity dispersion re- 
quired to make the BH more luminous than its central halo 
and also give a significant boost over the background, i.e. 
cq. I I25I I. The upper three dashed-dotted (blue) curves give 
the velocity dispersion of the parent halo, expression II26I I, 
from which the BHs have been ejected. Clearly a very small 
velocity dispersion is required which is more likely to be 
satisfied at higher redshifts. The bottom lower dashed (red) 
curve shows the evolution with redshift of the primordial 
dark matter velocity dispersion. 



observable diffused background flux. The total flux 
is given by the integral 



M z 



L{M, z) 
47rr(z)^ 



N{M,z)dMdV{z) 



(27) 



where M can be either the BH mass (Mbh) or the 
halo mass (Mhaio), r{z) = Rh — ^Ij^^ with 
Hubble radius Rh = 4000Mpc, and N{M, z) is the 
number density of the BHs [or haloes in the calcu- 
lation for NFW haloes] and the luminosity of a sin- 
gle BH L{M, z) (or the parent NFW haloes at z) is 
given by ((2T|) [or ((24l) ] and the volume element is 
dV = sinip r{z)'^dr(z)dip dip. 

The physical number density of the BHs is as- 
sumed to follow the Press-Schechter formalism 
(Press & Schechtcr 1974, Bower 1991), 
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halo 
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(28) 
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where n is the power spectrum index — 2 < n < — 1. 
For BHs, the above expression (|28p has to be mul- 
tiphcd by the relative time a BH spends at apapsis, 
i.e. A<df/^o where is the age of the Universe. 

The time spent at apapsis, (fT7|) . is itself a func- 
tion of BH mass. Thus only the fraction At up /to of 
the ejected BHs can be considered to be actually on 
their apapses passage in a Hubble time. We can eval- 
uate the total flux by performing the integrals in (1271) 
[multiplied by AtDp/io for BHs]. In using p7|) . we 
assume that the BH velocity Vbh is of the same order 
as the dark matter velocity dispersion (Tdm whose 
current value is 0.03 cm/s and increases linearly with 
(1 -I- z). We expect this assumption about the veloc- 
ity dispersion to be more valid at very high rcdshifts. 
For spectral index n ~ —1 and ^ 10^^ Mq and 
for the ensemble of BHs at their apapses passages 
and their central haloes (assuming NFW profiles) 
we obtain <i>NFW ~ 10^^ 7 cm^^sr^^ . The flux 
from the BHs is lower than this value, since we 
have evaluated our parameter i?cutoff by requiring 
Lnfw = ^BH- The BHs only spend a fraction of time 
at the apapsis which yields approximately $bh ^ 
10~^^ 7 cm~^sr~^.0 The flux would be atten- 
uated due to interaction of photons which however 
would affect approximately equally $bh and $nfw • 

We have assumed that only BHs produced in 3(t 
peaks of the density perturbation can undergo effec- 
tive mergers. We have assumed that all ejected BHs 
orbit their central haloes outside the virial radius; 
however, were this not the case, we do not expect 
any significant overall decrease in the flux which is 
already underestimated by our moderate choices of 
parameters and also by assuming that there is only 
one apapsis passage for a BH. We have also ignored 
the effect of multiple density-enhancement for a BH 
which is on its inward journey through an already 
high-density wake. 

6. Conclusion 

Black holes can be ejected from their host haloes 
due to anisotropic emission of gravitational waves 
in the merger of their progenitors. If ejected with 
velocities below the escape velocity they move on 
bound orbits around their host haloes. Since they are 
ejected from the centres of the haloes they move on 
radial orbits and their velocities come to zero at the 
apapses passages. Around their apapses, these BHs 

1 This corrects the typo in Mohayaco, Colin & Silk 2008 
where this ^bh was mistakingly written as ^nkw- 



have low velocities and move in a cold background 
and very high density wakes can form around them. 
If dark matter was to consist of self-annihilating par- 
ticles, these BHs would be powerful sources of high- 
energy 7-rays, both individually as resolved sources 
and collectively as diffuse background. The results 
here indicate that the globular clusters in the out- 
skirt of our halo or field galaxies in our local Uni- 
verse devoid of central BHs can have orbiting BHs 
which during their apapsis passages would produce 
flashes of high-energy 7-rays, although this effect 
is expected to be most signiflcant for massive ob- 
jects at high redshifts. The validity of dynamical 
friction formulae has been very rarely studied for 
radial orbits (Gualandris & Merritt 2007). The fact 
that there is no mass loss makes BHs a rare case 
for dynamical friction theory. Throughout this work 
we have assumed a homogeneous background and a 
constant-velocity approximation, both of these as- 
sumptions are questionable for the problem consid- 
ered here. The validity of these assumptions remains 
to be checked in high-resolution numerical simula- 
tions. 

We are grateful to Michel Henon, Joe Silk & David Wein- 
berg for contributions. We thank French ANR OTARIE for 
travel grants. 
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